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Abstract 

The discrete Chebyshev polynomials t n (x,N) are orthogonal with respect to 
a distribution, which is a step function with jumps one unit at the points x = 
0, 1, • • • , N — 1, N being a fixed positive integer. By using a double integral rep- 
resentation, we have recently obtained asymptotic expansions for t n (aN,N + 1) in 
the double scaling limit, namely, N — > oo and n/N — > b, where b £ (0,1) and 
a G (-00,00); see [Studies in Appl. Math. 128 (2012), 337-384]. In the present 
paper, we continue to investigate the behaviour of these polynomials when the pa- 
rameter b approaches the endpoints of the interval (0, 1). While the case b — > 1 is 
relatively simple (since it is very much like the case when b is fixed), the case b — > 
is quite complicated. The discussion of the latter case is divided into several sub- 
cases, depending on the quantities n, x and xN/n 2 , and different special functions 
have been used as approximants, including Airy, Bessel and Kummer functions. 
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1 INTRODUCTION 



The discrete Chebyshev polynomials t n (x, N) can be defined as a special case of the Hahn 
polynomials [U p. 174] 

Q n (x; a, (3, N) = 3 F 2 (-n, —x, n + a + /3 + 1; —N, a + 1; 1) 
(1.1) = ^ (-n) k (-x) k (n + a + P + l) k 

With a = (3 = in ([III]), we have 

t n (x, N) = (-l) n (N - n) n Q n {x; 0,0, N- 1) 

(I- 2 ) -(_-\\n(N- \ (-n)k(-x) k (n + l) k 

~ [ j ( )n ^ (-N+l) h klk\ ' 

see [H p. 176]. 

In a recent paper [5], we have studied the asymptotic behaviour of t n (x, N + 1) as n, 
iV — > oo in such a way that the ratios 

(1.3) a = x/N and b = n/N, 
satisfy the inequalities 

(1.4) - oo < a < - and < b < 1. 

In view of the symmetry relation [SJ eq.(8)] 

(1.5) * n (ar, N + l) = (-l) n t n (N -x,N + l), 

our study actually covers the entire real-axis — oo < x < oo or, equivalently, the entire 
parameter range — oo < a < oo. For a discussion of the asymptotic behaviour of the Hahn 
polynomials, see [1]. 



In the present paper, we shall investigate the behaviour of the polynomials t n (x, N + l), 
when the parameter b given in (11 .3p either tends to or tends to 1. The case b — > 1 turns 
out to be rather simple, since the ultimate expansions and their derivations are similar to 
those in the case when b is fixed; see (12.211) for < a < | and ("12 . 28[) for a < 0. 

To derive asymptotic approximation for t n (x,N + 1) when b — > 0, we divide our dis- 
cussion into several cases, depending on the quantities n, x and xN/n 2 . A summary of 
our findings is given in Table [1] below. Some explanation is needed for this table. By 
"small", "fixed" and "large", we mean, respectively, x = o(l), x G [S,M] and x ^> 0(1) 
as iV — > oo, where 5 is a small positive number and M is a large positive number. By 
"Airy", we mean an asymptotic expansion whose associated approximants are the Airy 
functions Ai(z) and Bi(z). In the same manner, by "Bessel" and "Kummer", we mean 
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Table 1: Asymptotic expansions in different cases when £ — > and N — > oo 



x 


small 


fixed 


large 


large 


Kummer 


Airy 


Bessel 


fixed 


Kummer 


Airy 


series (11. 21) 


small 


Kummer / series (11. 2 p 


Airy 


series (11. 21) 



asymptotic expansions whose associated approximants are, respectively, the Bessel func- 
tion J v (z) and the Kummer function M(a, c, z). The Kummer function used in this paper 
is defined by 

1 1 °° ( 1 s 

(i.6) %^) = fR M( '' v) = fiE^ 

which was introduced in [6, p. 255]. The case involving Bessel functions covers an earlier 
result of Sharapodinov [10J, where, instead of Bessel functions, Jacobi polynomials are 
used as an approximant. In the case when xN/n 2 is large and x is either fixed or small, 
the series in ( II. 2p is itself an asymptotic expansion (in the generalized sense [TTJ p. 10]) as 
N — > oo. Hence, there is no need to seek for another asymptotic representation. 

The arrangement of the present paper is as follows. In Section 2, we recall the major 
results in [8j; to facilitate our presentation for later sections, we also include here brief 
derivations of the expansions given in [8]. In Section 3, we consider the case when xN/n 2 
is small and show that the expansions are of Kummer-type. In Section 4, we consider the 
case when the quantity xN/n 2 is fixed; that is, when xN/n 2 is bounded away from zero 
and infinity. There are three subcases, depending on x being large, fixed or small. In all 
three subcases, we will show that the expansions are of Airy-type; see Table [TJ The case 
when both quantities x and xN/n 2 are large is dealt with in Section 5, where we show 
that the expansion of t n (x, N + 1) can be expressed in terms of Bessel functions. The 
final section is devoted to two remaining cases, namely, (i) when xN/n 2 is large and x is 
either fixed or small, and (ii) x < 0. 



2 RESULTS IN REFERENCE E 



In [8], we have divided our discussion into two cases: (i) < a < ~, and (ii) a < 0. As 
mentioned in Section 1, by virtue of the symmetry relation (II. 5p . these two cases cover 
the entire range — oo < a < oo. 

In case (i), we started with the integral representation 

(2.1) t„ (x, N + i) = <=£ ^ + 2> +1 , f I - V*«^d t , 

2m T(n + l)T(N - n + 1) J — 1 

where 

(2.2) f(t,w) = 61n(l - t) + (1 - b) hit + alnw - aln(w - 1) + 61n[l - (1 - t)w), 
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and the curve 71 starts at w — 0, runs along the lower edge of the positive real line 
towards w — 1, encircles the point w = 1 in the counterclockwise direction and returns to 
the origin along the upper edge of the positive real line. As a function of two variables, 
the partial derivatives df/dt and df/dw vanish at (t,w) = (t±,w±), where 

, , 2 - 2a - b 2 + bVb 2 - 4a + Aa 2 b + Vb 2 - Aa + 4a 2 

( 2 - 3 ) U = S7i \/i , u\ ' w + 



2(l-a)(l + 6) ' ^ 26 



and 



, , 2 - 2a - b 2 - b^b 2 - 4a + 4a 2 6 - V6 2 - 4a + 4a 2 

2.4) t_ = ; — , W_ = . 

v ; 2(1 -a) (1 + 6) 26 

For each fixed w 6 71, we now find a steepest descent path of the phase function f(t, w) 
in the variable t, which passes through a saddle point t Q (w) depending on w. (Note that 
not only the saddle point t (w), but all points t on this steepest descent path depend on 
w.) The function f(t (w),w) is a function of w alone. To find the relevant saddle point 
to(w), we solve the equation df(t,w)/dt = 0, and obtain 

2w-l + ^/l + Ab 2 (w-l)w 
2(1 + b)w 

It can be shown that 

(2.6) to(w+) = t+, to(w-) = t- 
Define the standard transformation t — > r by 

(2.7) T 2 = f(t (w),w)-f(t,w); 

see p!Tl p. 88]. Note that we have r = —00 when t = 0, and r = +00 when t = 1. 
Furthermore, this mapping takes t = t (w) to r = 0. Coupling (12.11) and (12 .7p gives 

*.(»," + !) r (" + Ar + 2 ) 



2« r(n+l)r(AT-n + l) 

^ ' /•+°° /" 1 rl + 

x / / -^—e N ^ w ^e- NT —dwdr, 
J -00 dr 

where 71 is the integration path of the variable w in (12.11) . Note that the first variable of 
/ in (I2.8P is now t (w), instead of t. Hence, the phase function is a function of w alone. 
Setting df(t (w),w)/dw = 0, we obtain the saddle points 



6 ± Vb 2 - 4a + 4a 2 

(2.9) w ± = ; 

cf. (l2.3p and (12. 4p . Motivated by (12. 2p . define the new phase function 

(2.10) ip{u) = a In m — a ln(w — 1) + rju. 
The saddle points of if) are given by 



(2.11) u ± = . 
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To reduce the double integral in (12. ip into a canonical form, we define the second mapping 
w — y u by 



(2.12) 

with 
(2.13) 



f(t (w),w) = i/j(u) +7 

= alnw — aln(u — 1) + rju + 7 



u(w + ) = U. 



u{wJ) = u + , 



where 77 and 7 are real numbers depending on the parameters a and b in (I2.2p . From 
flZ32D , we have 



(2.14) 



dw dip /df(t (w),w) 
du du / dw 



r)(u — u+)(u — U-)w(l — w) [(2a — 1) — y/T+ Ab 2 w 2 — Ab 2 w 



-2b 2 u[u — l)(w — w + )(w — W-) 



for u 7^ u±, and by L'Hopital's rule, 



(2.15) 



dw 
du 



^rf + 4arj(l - 2a) (1 - a)(-r)) 
- 4a + 4a 2 



1/2 



for u = u±. With the change of variable to — > w defined in (12.121) . the representation in 
(12. 8 p becomes 



* n (z, JV + 1) 



T(n + N + 2) 



(2.16) 

where 
(2.17) 



2m Y{n + l)r(iV - n + 1 
"+°° /■ a(u,r^ 



3 7V 7 



X 



00 1/7. 



u-1 



N^u)-NT^ dud 



h(u, t) 



u — 1 dto dt 
w — 1 du dr 



and 7„ is the steepest descent path of ip(u) in the u-plane. An asymptotic expansion, hold- 
ing uniformly for a G [0, |], was then derived in [8] by an integration-by-parts technique. 
To state the result, we define recursively h (u,r) = h(u,r), 



(2.18) 
and 



(2.19) h l+1 (u,r) = (u-l) 

Furthermore, we write 
(2.20) 



hi(u,r) = ai(r) + bi(r)u - (u - u-)(u - u + )gi(u,r) 

dgi(u,r 



gi(u,T) + u- 



du 



1 = 0,1,2,.... 



j=Q j=0 
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The resulting expansion takes the form 

t . n \ AT™ ( OO 



-l) n T(n + N + 2)e N i f , ^ 

— K — ; ^_ <^ M(aN + 1, 1, r/N) V 

+ 1 )T(N-n + l )Vn\ K U 



tn ^ N+l) Vill , _ „ .... | ' ' ""f^(//A7 

(2.21) ' ' 

+ M'(aN + l,l, V N)J2-^ )n , 

where M(a, c, z) is the Kummer function defined in f ll.6p . and the coefficients q and di 
are given explicitly by 

- 1 - 1 

(2.22) q = ai- m ,2mT(m + -)r] m , d t = ^ bi- m ,2mT(m + -)r] m . 

m=0 m=0 

Case (ii) was dealt with in a similar manner. We started with the double-integral 
representation 

( , 23) t n {x,N + i) = H>I I ( "^ + 2) +n f f - V^w, 

27T« r(n + l)r(iv - n + 1) Jo J l2 w 

where 

(2.24) /(*,«;) = Mn(l - t) + (1 - b) hit + aln(-w) - aln(l - w) + Mn[l - (1 - t)w], 

and the integration path 72 starts at w = 1, traverses along the upper edge of the positive 
real line towards w — 0, encircles the origin in the counterclockwise direction and returns 
to w — 1 along the lower edge of the positive real line. 

Now we follow the same argument as given in Case (i), and use the same mapping 
t — > r defined in (12. 7p . except with f(t,w) replaced by f(t,w). The result is 

t (xN + 1) J~ 1)n r ^ + iV + 2 ) 

(2-25) +oc 

x / / l^M^i^. 

Because of the shape of the contour 72, it turns out that the phase function f(t (w),w) 
has only one relevant saddle point, namely see (12.91) . As a consequence, we define the 
mapping w — > u by 

(2.26) f(t (w),w) = aln(-ti) -it + 7 
with 

(2.27) = a, 

where 7 is a constant depending on the parameters a and 6 in (II. 3p . The final expansion 
is in terms of the gamma function, and we have 

00 



(2.28) t n (x, AT + 1) ~ T{n + mN _ n + m _ aN + 1) ]L ^ 

where q are constants that can be given recursively. 
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3 KUMMER-TYPE EXPANSION 



As b 0, some of the steps in Section 2 are no longer valid. Let us first examine the 
mapping £ — > r given in ( 12. 7ft . In this mapping, we have used the fact that for each fixed 
w £ 7i, the saddle point £ = t (w) in (12 ,5p is bounded away from £ — and £ = 1, and 
the steepest descent path from t = to t = 1, passing through £o(w), is mapped onto 
(— oo, oo) in the r-plane. However, from (12.51) . we have for any fixed w 

(3.1) £ (» = 1 - b + O (b 2 ) as 6^0; 

that is, as b — > 0, to(w) approaches the branch point t = 1 in the t-plane. Thus, the 
mapping (12.71) is no longer suitable in this case. Since we are more interested in the 
neighbourhood of t = 1, where the term Mn(l — £) in the phase function (12.21) becomes 
singular, we introduce the mapping 

(3.2) f(t,w) = Mnr-r + A, 

where the constant A does not depend on t or r (but may depend on w)\ see (13. 4p 
below. To make the mapping t — » r defined in ( 13. 2p one-to-one and analytic, we prescribe 
t = to(w) to correspond to r = b, which is the saddle point of 61nr — r + A; i.e., 

(3.3) r(£oH) = b. 
This gives 

(3.4) A = f(t (w),w) -b\nb + b. 

Note that we have r = when t = 1, and r = +oo when t = 0. Furthermore, from ( 13. 2p . 
we have 



d£ b/r-1 (t -&)£(!-£) [1 - (l-t)w) 



dr df(t,w)/dt T(l + b)w(t-t+(w))(t-t (w))' 

dt f4H(i-4H)[i-(i-4WH^ 12 



r^b, 



dr [ by/1 + Ab 2 (w - l)w 

where we have used L'Hopital's rule for r = b and 

+ / , 2w - 1 ± JT+ Ab 2 (w - l)w 

to {w) = 27TT% ; 

cf.flZSD. Coupling (j2DJ and (E2D, we have, instead of (ESD, 

(-l) n+1 r(n + + 2)e n &-" 
£ n (x,iV+l)- — r(n + 1)r(A r_ n + 1) 

(3-6) , +00 



d£ 



. e JV/(*0 (t«),t«) e A r (b In r-r) _^ dwdr _ 



71 w — 1 dr 

Next, let us examine the mapping (12.121) . We note that this mapping still works in 
the present case; the only difference is that for fixed n/N G (0, 1), the constant 7] in the 
mapping (12. 12[) is bounded away from 0, whereas when n/N — y 0, rj approaches 0. More 
precisely, we have rj ~ —b 2 . This can roughly be seen from the saddle points (I2.9p and 
(12. lip , together with their correspondence relation (12.131) under the mapping w — > u given 
in (I2.12p . To prove this rigorously, we give the following lemma. 
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Lemma 1. Let rj be the constant defined in the mapping A2.12\) . If a/b 2 = o(l) (or, 
equivalently, xN/n 2 = o{l)), then we have 

(3.7) v = -b 2 + o(b 2 ) asb^O. 
Proof. From (12.51) . we have for any fixed w 

(3.8) U-) = T ^-%^-^^ + 0( () V(--l) 3 ) as^O, 

where the term O (b 4 w(w — l) 2 ) holds uniformly when either w or 1 — w is small. If 
a/b 2 — > 0, from (12. 9p we have = 1 — w + —> + and both w± bounded. Thus, 

to(wJ) = l-b + b 2 {l - b)w- + O (fcW ) , 

(3.9) 1 h 2 

t o(w+) = ^ - + O (&V) . 

To obtain 77, we use (1232]) and (1213]) . First, substituting O in (Q gives 

1 — -in 

(3.10) f(jt (w+),w+) - f(jt (w-),w-) = 2a In - b 2 + O (6 3 ) . 

Here, we have made use of the fact that if w± lie on the upper edge of the cut along the 
interval < w < 1, then we have \n(w± — 1) = ln(l — w±) + iri. Similarly, if w± lie on 
the lower edge of the cut, then ln(«;± — 1) = ln(l — w±) — ni. Next, we have from (I2.10p 
and (1211]) 



, 1 —u+ / 4a 

3.11 ${u-) -Mu+) = 2aln - + + — . 

u+ V V 

Note that u + + U- = 1 and rj < 0. Formula (13.71) now follows from a combination of 
( 12121) . (12~T3|) . ( 13101) and (131TD . □ 

What Lemma 1 says is that if a/b 2 — > and x = 0(1) or i > ^(1), then we have 
\r)\N ~ o 2 iV = xb 2 /a — > oo, i.e., r/A^ is large. 

Coupling (I3.6p and (I2.12p . we have 

r_i)n+ir(n + N + 2)e n b~ n e N ~< 



t n (x,N + l) 



(3.12) 



where 



2mT(n + l)r(iV - n + 1) 



x 



71 



!}M e ^)^>»r-r) dll , dT 
U — 1 



/o i o\ \ u-\ dw dt 

3.13 h{u,T) = 

w — 1 dw dr 

dw/dw is given by (12. 14j) and (12. 15[) and dt/dr is given by (13. 5p . Following the same 
integration-by-parts procedure outlined in Section 2, we let h (u,r) = h(u,r) in (I3.13p . 
and define hi(u,r) recursively by (I218p . (12.191) and (12.201) . The final result is 



8 



t n (x,N + l) 
(3.14) 

where 
(3.15) 

(3.16) 

and 
(3.17) 



(-l) n+1 r(n + N + 2)e n b~ n e N ^ 



x 



T(iV - n + 1)^"+! 

p-1 



M(aiV + 1, 1, V N) £ j^y + M\aN + 1, 1, rjN) ]T ^ + £ 



A r 



n+l /-co 



r(n + 1 



i \ n -Nt A T(n + m + l) 

a; r r e dr ~ > a/ m — : : . 

iy 1 ^ l ' m T(n + l)N m 

oo 



iuri+1 poo ^ 
v ; JU m=0 



r(n + m + 1) 
r(n + l)A m 



A r 



n+l 



27rzr(n + i)(^iv)p y 



^Vr) e JV^(«) e JV(6 In T-r) dMdr 
, " 1 



Note that since 6 = n/N — > 0, the series in ( 13. 15[) and ( 13 . 1 6[) are asymptotic. More- 
over, since r)N is large in this case, (13. 14j) is indeed a compound asymptotic expansion as 
N ^ oo. 



Now, we consider the subcase in which a/b 2 — > and x = o(l). In this case, the 
expansion in ( 13. 14|) is still asymptotic as long as \rj\N ~ & 2 A/" — > oo or, equivalently, 
n 3> 0(y/~N). If n = 0(y/N), then the series in ( 11 .21) is itself an asymptotic expansion 
as — )• oo. This completes our discussion of all three cases listed under the condition 
u xN/n 2 small" in Table 1. 



4 AIRY- TYPE EXPANSION 

For the case a/b 2 £ [7, M], where 7 is a small positive number and M is a large positive 
number, the saddle points w± in ( 12 .9f) are bounded away from the singularities w = 0, 
w = 1 and w = 00, and coalesce with each other when a/6 2 — > ^yz^y. Therefore, to derive 
an asymptotic expansion uniformly for a/b 2 £ [7, M], we need the cubic transformation 
(see (14. 6p below) introduced by Chester, Friedman and Ursell [2] . The resulting expansion 
is in terms of the Airy function Ai(-). 

Following the same argument as in Section 3, we again use the mapping from t — > r in 
( 13. 2p . and start with the integral representation ( 13. 6 j) 



2ni r(n + l)r(JV - n + 1) 

l 4 - i J /-+oo 1 J f 

X / / _J_ e Nm( W ), W ) e N(Mn T - T )M d 



j 7l w - 1 dr 



where the integration path is described in the line below (12. 2p . To proceed further, we 
divide 71 into two parts, and denote the part in the upper half of the plane by 7^ + \ and 
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the other part in the lower half of the plane by 7^ >. Recall from f 12 . 2 jl that the phase 
function in (14. ip is given by 

f{t {w),w) =61n(l - t (w)) + (1 - h) \nt Q (w) + alnw 
- a\n(w - 1) + b\n[l - (1 - t (w))w], 

which has a cut (—00, 1] in w-plane. Put 

J{t {w),w) =61n(l - *oH) + (1 - 6) IntoH + alnw 
1 ' J - aln(l - w) + 61n[l - (1 - * o («0H> 

which has cuts (— 00, 0] and [l,oo). From (14. 2p and (14. 3p . we have 

(4.4) f(t (w),w)=J(t (w),w) + airi 
for w G 7 ( - + \ and 

(4.5) f(t (w),w) = J(t (w),w) - ani 
for w G 7 ( -~* ) . We now make the standard transformation 

(4.6) J(t (w),w) = ^u 3 -(u + A J 

with the correspondence between the critical points of the two sides prescribed by 

(4.7) w+ y/l, w^^-^/C. 

If C > 0, then ±\/C are both real; if ( < 0, then ±vT are complex conjugates and purely 
imaginary. The values of A and ( can be obtained by using (14.61) and (14 .7p . We also have 

dw (u-y/£)(u + V£)w(l - w) \(2a - 1) - VI + 46 2 w 2 - 4b 2 w\ , r- 

u t ±VC, 



du — 2b 2 (w — w+)(w — w. 

dw f2v^o(l-2a)(l- 



(4.8) . 

dw [ oV^ 2 - 4a + 4a 2 

where we have again used L'Hdpital's rule for it = ±a/C Let us first consider the case 
C < 0, and deform the image of the contour 7W under the mapping w — >■ u defined in 
(14. 6p to the steepest descent path of |m 3 — (u + A in the u-plane which passes through 
\fQ. We denote the path by C<i- Similarly, we deform the image of the contour and 
denote the steepest descent path passing through — \/C by C3; see Figure 1. Next, we 
consider the case ( > 0. Note that here the saddle points ±a/C are real, and that the 
contours C2 and C3 pass through both of them; see Figure 2. Clearly, in both cases, C3 
is the reflection of C2 with respect to the real axis in the w-plane. 

Let C\ denote the dotted curve shown in Figures 1 and 2, and recall the identities 
e iaNn _ qos aiW + ismaN7i and e~ mNn = cosaNir — isinaNrr. A combination of 04. ip . 
( fOll . f l4~5l) and flO} then gives 

, Ar , (-l) n T(n + N + 2)e n b~ n e NA 

t n ( X , N+l) = V ' \.s v(M -Ts 

I (n + lji (iV — n + 1) 



(4.9) x f °° [ h (u,T)e N (^)e N ^-r) dudT 

2tt« Jq J Cl 

sin aNn r+ °° 



2tt 



h (u,r)e N ^ u ^ u )e N ^ T -^dudr 
ii </c 3 -c 2 
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where 

,,.„, , , s 1 dt dw 

4.10 h (u,r) = r^-p: 

w — 1 dr au 

dt/dr, dw/du are given respectively by (13. 5 j) and ( 14. 8p . 




4.1 When £ is bounded 

Following the standard integration-by-parts procedure [TTJ p. 368], we define recursively 

(4.11) = cli(t) +bi(r)u- (u 2 -()gi(u,T), 

(4.12) hi +1 {u,r) = — — — . 
Furthermore, expand 

oo oo 

(4.13) a,(r) = J>y^, &,( r ) = ^6,^'. 

j=Q j=0 

From gS}, pi) . (Q2]) and (Qgjl . it follows 
. (-l) n r(n + iV + 2)e n 6-"e JVA 

f " (l - iv + 1 > = r ( AT- n + i)U' 



(4.14) 



( Ai(iV^C) £ ^ - Ai'(iV^C) £ 
+ sin aN-n [ Bi(iV 2/3 C) ^+1/3 ~~ 

Bi'(iV 2 / 3 C) N l+2/3 ) + 6 P 
1=0 1=0 



where 



(4.15) 



£ 



P 



h r {u,T)e N (^'- i ")e r "- b '"-^AuAT, 



JCx 



2iriNPT(n + 1) 
sin (aN7r)N n+1 



27rziVPr(n + l)7 7 C3 _ C2 
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and the coefficients are given by 

^ n+1 f° ^ n -Nrs ST* T(n + m + l) 

Ci = —, r / ai(T)T e dr ~ > ai m —— — - — , 

r(n + l )Jo lK } ^Tn + liV™' 

(4.16) m -° 

Nn+1 ~ r(n + m + l) 

di = —, r / bi{T)T n e dr ~ > km—, r -. 

r(n + l)7 iy 1 ^ '.«r(n+l)JV* 

In the present case, a/b 2 G [5, M] or, equivalently, w± are bounded. Thus ( is also bounded 
in view of (14. 7ft . Therefore, it is easy to prove that (14. 14ft is asymptotic; for details, see 
HD P-371-372]. 

4.2 When ( — > — oo 

When a/b 2 — > oo, by (12. 9 p the saddle points w± approach | ± zoo, respectively, along the 
line Re w = |. In view of (14. 7p . this is equivalent to saying that ( — >• — oo. To prove that 
(I4.14p is also asymptotic in this case, we rewrite the expansion in (I4.14p as 

, (-l) n T(n + N + 2)e n b- n e NA 
+ r ( iV-n + l)U 

\ *r [Ai(N 2 / 3 C) P ^ ci Ai'(iV 2 / 3 C)^ d x 
4 - 17 x cos aNn I \ —L— \ ±± y —L 

+ SmaiV7r {-l^lsJaJty- N 2 / 3 l.JaW l ) + f 
where e p is the same as in (I4.15P and 
(4.18) q = a l ci and di = a l d[. 

Thus, it is sufficient to first prove the boundedness of the coefficients q and di, and then 
establish the asymptotic nature of the error term in (I4.17p . namely, to prove that there 
exist positive constants A p , A' , B p and B' p such that 



gl . < K Ai(N 2 / 3 () | A' p Ai(N 2 / 3 () 



(4.19) 



where 



(4.20) 



v 



{aN)P N 1 / 3 (aN)P N 2 / 3 



e 



, 2| j?„ Bi(iV 2 / 3 C) g Bi (jV 2 / 3 Q 
|£p| -(aAT)P jVV3 + (aiV)P iV 2 / 3 

x cos(aNir)N n+1 r+ °° 



p 2niNPT(n + 1) J J Cl p 



s 2 



sm(aN-K)N n+1 f +OD 



> 2niNPT(n + l)J Q J Cs _ C2 p 



■ > 1 /2 

and M(z) = Bi(z) = [Ai 2 (z) + Bi 2 (,2)] for z < 0. Note that in the present case, ( < 0. 



h p (u,r)e N (^ u )e N ^-^dudr, 
h p (u,r)e N ^ ui - iu )e N ^ bhlT -^dudT, 
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From (14. 8p . we recall that dw/du depends on (. The value of ( is obtained by solving 
the two equations gotten from (14.61) with w and u replaced respectively by w± and ±\/C; 
see [HI p. 367]. Thus, ( and hence dw/du both depend on the parameters a and b in (11.31) . 
As a and b approach zero, ( may tend to infinity. When the saddle points u± = i-^/C 
are bounded, the mapping w — >■ u defined by the cubic transformation (14. 6p is analytic, 
and its derivative dw/du is bounded. However, in the present case, the saddle points 
u± = ±\/— (i go to infinity as ( — > — oo. Hence, the coefficients aj(r), 6/(r) and the 
function hi(u,r) given recursively in (14.111) and ( 14. 12|) may all blow up as ( — > — oo, 
where h (u, r) is defined by (14.101) . Since the coefficients q and di in (14. 14j) are related to 
aiij) and hij) via (14.161) and (14.131) . to prove that the expansion in (I4.17P is asymptotic 
when C — oo, we must first give estimates for the coefficient functions a;(r) and bi(r). 
To this end, we shall adopt a method introduced by Olde Daalhuis and Temme [5J. To 
begin with, we define 

(4.21) Ao (u,0 = -^-, B (u,Q 



U z — Q U z — Q 

Using (14. lip and the Cauchy residue theorem, it can be verified that 

a o( T ) =77— / ho(u,r)A {u,C)du, 



r 



(4.22) 2 f 

b o( r ) =77— / ^(^^^(^Odw, 

ZTTl 



r 



where T is the contour consisting of two circles, centering at both with radius R, 

where R could be as large as possible until the circles reach the singularities of h Q (u,r) 
in the w-plane. 



We note that w = w± are removable singularities of dw/du. However, dw/du blows 
up at the points w = w±e 2km , where k 7^ is any integer; see (14. 8p . To find their image 
points in the w-plane under the mapping (14. 6p . we take Wk = w + e 2km as an example and 
denote its image point by Uk- The image point of W-e 2klTl can be treated in a similar 
manner. By (14. 6p . we have 

(4.23) 7(to(w k ),w k ) = -ut-Cu k + A. 
From (14. 6 p and (14. 7p . we also have 

(4.24) J(t (w + ),w+) = ^u 3 + - (u + + A. 
Subtracting (I4.24[) from (I4.23P gives 

(4,5) 2^<-<p(44-4 

Since < a < \ and -( > 0(1), from f l4T25|) it follows that \z k \ ~ V^a 1 / 2 ^)- 3 ^. 
Therefore, there exists a constant < Cq < 1, independent of a, b and (, such that the 
interior of the two circles with centres at u = ±y/—£i and radius R = CQ^Ja{—C ! )~ l ^ is 
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free of the singularities of ho(u, r) in the u-plane. Since ho(u, r) is now analytic inside the 
contour T, there exists a constant Ch such that 

(4.26) \h (u,r)\ < c h h{r) 

for u in the domain enclosed by the contour T, where h(r) denotes the maximum of the 
two functions |/io(±v^~C^ r ) | • Note that as functions of r, ho(±y/—(i, r) are analytic in 
the neighbourhood of steepest descent path in the r-plane. 

We further introduce rational functions A k and B k , k = 0,1, 2..., defined recursively by 

A k+1 (u,C) =-^—-^-A k (uX), 

(4.27) 1 d^ 

B k +i(u,() =— — --7-B k (u,C), 
u — Q (xu 

where Aq and Bq are given in (I4.2ip . By induction, we can show that A k and B k are 
expressible as 



[(fc+l)/2] „1 „,fc+l-2 

A k {u,()= 

(4.28) 



„2 „,*-« 



£ fc (u,C)=^ 



t=0 



U 2 - C) 



where c ki and c 2 ,^ are constants independent of u and (. As in (I4.22p . by Cauchy's 
theorem, we have from equations (14. lip and (I4.12p 

a k {r) =— / /i fc (M,r)A (M,C)dM 



r 



hk-i(u,T)Ax(u,()du - — I (afc-i(r) + b k -x{r)u) A x (u, Qdu 



2lli ./ r 27T2 j r 



( 4 - 29 ) =^ / /*_i(u,r)Ai(u,C)du 



1 



2tu 



/i (u,r)v4 fc (u,C)dM, 



r 



where we have used integration-by-parts to derive the second equality. The second term in 
the second equality vanishes because (a k -i(r) + b k -\{r)u) A\(u, () is O (u~ 2 ) as \u\ — > oo 
and all poles of that function lie inside T; see (14.28)) . Similarly, we also have 

(4.30) b k (r) = -L / h (u,r)B k (u,C)du. 
Using (14.281) . it is easy to obtain the estimates 

(4.31) \A k (u,C)\ < C k a- k -^ 2 {-0 1/A and \B k {u,()\ < C k a~ k , 
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for u on and inside the contour Y and — £ — > oo. Here and thereafter, C k is used as a 
generic symbol for constants independent of u, (, a and b. Substituting (I4.3ip into ( I4.29[) 
and f)4.30p gives 

(4.32) \a k (r)\<C k h(r)a- k , 
and 

(4.33) |6 fc (r)| < C k h(r)a- k+1 / 2 (-Cr 1/4 . 

Therefore, a l ai(r) and a l bi(r) are both bounded for r in the neighbourhood of steepest 
descent path in the r-plane, and of course for r in the neighbourhood of r = 0. Thus we 
have the boundedness of the coefficients q and d\. 

To estimate h k (u,r), we use the rational functions R k (u,w,(), k = 0,1,2..., defined 
recursively by 



(4.34) »"•" d 

R k+1 (u,w,() =— 2 — -—R k (u,w,()- 
u z — Qau 

These functions were also introduced by Olde Daalhuis and Temme [5]. They showed by 
induction that R k (u,w,() can be written as 

k — 1 min{i,fc— l—i} „ j_j 

(4.35) R k (u,w,0 = J2 E 7 2 Tvw ' k = l,2..., 

where Cy do not depend on u, w and C- Similar to f)4.29p . we have 

(4.36) h k (w,r) = — [ h (u,r)R k (u,w,()du, 

where T is the same contour used in (14.291) and w lies inside two disks centered at ±V — 
and with radius |co\/a(— C) -1 ^ 4 - It is easy to verify from (I4.35j) that 



(4.37) I^K^OI^^a-^/^-C) 1 / 4 , 
and from (14.361) that 

(4.38) \h k (w,r)\<C k a- k h(r). 

Substituting (14T38I) into (14T201 gives (14TT91 : for details, see [5j p.311-312]. Note that to 
make the expansion (14. 17j) asymptotic, we require x = aN to be large. 



5 BESSEL-TYPE EXPANSION 

In the case of Hahn Polynomials Q n (x; a, (3, N) given in (jl.ip . Sharapodinov [10] has 
given an asymptotic formula involving Jacobi polynomials when the parameters satisfy 
a, (3 > \ and 2 < n < cy/N, where c is a positive constant. The values of the variable 
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x are also required to be large; more precisely, x > eN and e > is a small number. 
Although discrete Chebyshev polynomials t n (x,N) given in (jl.2p is a special case of the 
Hahn polynomials, the values of the parameters are a = (3 = 0; that is, Sharapodinov's 
result does not include our case. However, since the leading term in the uniform asymp- 
totic expansion of the Jacobi polynomials is a Bessel function (see [TJ p.451]), the work of 
Sharapodinov did inspire us to look for an asymptotic expansion for t n (x, N + 1) involving 
Bessel functions, when the parameters a and 6 in (I1.3P satisfy a/6 2 — > oo and the variable 
x is large. Our method differs completely from that of Sharapodinov. 

Returning to (12. ip . and making the change of variable v = 1/w, we have 

(b .d tnix , N+1) = (zg ?".;% +2 > r I -i-™,*, 

2tu T(n + 1)T(N - n + 1) J Q J 73 v(l-v) 

where 

(5.2) f(t, v) = 61n(l - t) + (1 - b) hit - aln(l - v) - blnv + b\n[t + v - 1]; 

the curve 73 starts at t> = +00, runs along the lower edge of the positive real line towards 
v = 1, encircles the point v = 1 in the clockwise direction and returns to +00 along the 
upper edge of the positive real line. Note that since e Nblnv = v n , there is no need to have 
a cut from the origin to infinity. 

The saddle point of f(t,v) in the t-plane is given by 



(5.3) t {v) = 



v 2 



2(1 + 6) 

cf. fl2.5p . where the v -plane is cut along two line segments joining to the two con- 
jugate points 2b 2 ± 26^1 — b 2 i, and the branch of the square root is chosen so that 
a/46 2 — 4b 2 v + v 2 ~ v as v — > 00. From (15. 3p . we have 

t (v) = l + O(b) as \v I ~ cb, 
to(v) ~ 1 — 6 as f ^> 6, 

where c is a positive constant, and for i;<6, 

^ 5 5 ^ <o(v) ~ 1 - 2 U as w ^ + , 

t (f ) -^1-26 as v -)■ 0~, 

where + and 0~ mean limits approaching from the right-hand side of the cut and the 
left-hand side of the cut, respectively. Moreover, easy calculation shows that for any v, 
to(v) is not a real number on the cut (1, 00) in the t-plane. Following the same argument 
given prior to (13. 2p . we introduce the mapping 

(5.6) f(t,v) = b\nr -r + A 

with the correspondence between the saddle points t = t (v) and r = 6 given by 

(5.7) r(t (v)) = b. 



16 



Coupling (I5.6P and ( 15 .7p yields 



(5.8) A = f(t (v),v) -b\nb + b. 

The zeros of df(t,v)/dt = are given by 

2-v± - Ab 2 v + v 2 



(5.9) t*(v) 



2(1 + 6) 



where is the relavant saddle point to(t>) given in ( 15.31) . By straightforward calcula- 

tion, we have 



dt 6/r-l (r - b)t(l - t)(t + v - 1) 



(5.10) 



dr df (t, t(1 + &)(*- - to 0)) 

dt ft + W(l-to + M)(4W + ^-l)' 12 



dr [ b\/Ab 2 - Ab 2 v + f 2 



= 6. 



Note that we have r = when t = 1, and r = +oo when t = 0. For any fixed we 
can deform the original interval of integration < t < 1 into a steepest descent path r„, 
passing through t (v). Also note that t r (v) is not on T v , unless < v < 1. Moreover, 
if < f < 1, then < t (f ) < t^ (f) < 1 and r„ is the real interval [0, 1]. In our case, 
there is only one point on the path 73 in the w-plane (see ( 15. ip ). where it crosses the 
real line. Let us denote this point by vq- When v = vq, the mapping which we have 
introduced in (I5.6P becomes singular at the point t r (v), since dt/dr in (15.101) blows up. 
However for this particular case, we only need to slightly modify the path by replacing 
part of the original path near this point by a small half circle as shown in Figure The 







• • — •— > 



Figure 3: The indentation of the integration path in the t-plane for < v < 1. 

corresponding integration path in the r-plane following the mapping ( 15 .6p also needs to 
be modified. However, this small modification on the integration path will not affect the 
following argument and calculation. With this in mind, we will simply ignore this par- 
ticular case, and proceed with the assumption that the integration path in the r-plane is 
always [0, +00) for all w, and that dt/dr will not blow up in the neighbourhood of the path. 

Thus, from (15 .ip and (15.61) . we have 

-l) n+1 r(n + iV + 2)e n b~ n 



t n (x,N + l) 



2vri T(n + 1)T(N - n + 1) 

(5 - U) ,+oo , . df 



X 



J yi 



v(l — v) dr 
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Here, we rewrite the phase function f(to(v),v) as 

(5.12) f(t {v),v) = (l-b)lnt (v) + bla(l-t {v))-aln{l-v)-bln( 



t (v)+v-l' 



Recalling the statements following (15.21) and (15.31) . we know that there are only two cuts 
in the v -plane: one along the infinite interval [1, +oo) and the other along the bent line 
joining the conjugate points 2b 2 ± 2byl — b 2 i and passing through the origin. To find the 
saddle points of f(t (v),v), we set 

#oM =0) 

dv 

and obtain 



. . b 2 =F biy/Aa - Aa? - b 2 

(5 - 13) v± = ■ 

Since (15. ip is obtained from (12. ip by making the change of variable v = 1/w, (15.131) can 
also be derived from (I2.9p . Note that in this case, b 2 /a — > 0. Furthermore, since a < |, 
the quantity inside the square root is positive. Hence, v± are distinct, and approach v — 0. 

Define 

g(u) = m{u ), 

u 

where m > is some constant to be determined. The saddle points of g(u) are 

(5.14) u± = ±i. 
Make the transformation 

f(to(v),v) =g(u) + -f 

(5.15) 1 

= m(u ) + 7 

u 

with 

(5.16) u ( v +) — u -i u{vJ) = u + . 

Note the fact that to(v±) = t±, where t± are given in (I2.3P and (12.41) . This can be seen 
from (12. 6 p and the change of variable v = 1/w that we have made. Thus, substituting 
(IQBjl into (15TT5]) gives 

1, 1-6 1 , b 2 
7 = - In H — o In 



(5.17) ' 2 1 + b 2 1-b 2 

= b\nb-b + 0(b 3 ) 

and 



1 r 7N by/Aa - Aa 2 - b 2 by/Aa - 4a 2 - b 2 

m = — < ( 1 — b) arctan — a arctan — — 

2l v ; 2-2a-6 2 2a-2a 2 -6 2 



yj '4a — 4a 2 — b 2 

(5.18) -2b arctan- 

2 + b - 2a 

b arctan * / — - — , 
1 — a 
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when b 2 /a — > and b — > 0. Moreover, from (I5.15P we have 

dv m(u -i)(u + i)v(l - v) Ul - 2a)v + ^4b 2 - Ab 2 v + v 2 



(5.19) 



du 

dv 
du 



u=u± 



2a(l — a)u 2 (v — v + )(v — i>_) 
-2m(v T ) 2 (l - v T )(l - 2a) ) 1/2 
6V4a - 4a 2 - b 2 J ' 



u ^ u±, 



u = u±. 



Here we have made use of the equality a/46 2 — Ab 2 v T + = (1 — 2a)v T . 

Coupling (15. lip and (I5.15p . the integral representation of t n (x, N + 1) becomes 

(-l) n+1 T(n + N + 2)e n b~ n e N " / 

^n(-£> "I - 1) " 

(5.20) 



x 



2?ri r(n + l)r(iV - n + 1) 

" / ^ T ) e Jy[w(^l)] e JV(6lBr-r) dwlr> 

o h « 



where 
(5.21) 



/l(it, r) 



it dt dv 
f (1 — t>) dr dw ' 



and the contour 7 starts from —00, encircles the point u = in the counterclockwise 
direction and returns to —00. For Z = 0, 1, 2, • • •, we define recursively 



(5.22) 
(5.23) 



hi(u,r) = o,(r) + ^+(l + i] 9l (! tl r) I 

u \ u z J 

hi+i(u,r) = 



u d { qi(u,t) 



mdu \ u 



where ho(u, r) = h(u, r) given in ( 15.211) . Furthermore, we expand a/(r) and h(r) at r = 0, 
and write 



(5.24) 

It is easy to see that 
(5.25) 



a i ( r ) = ^2 ai >i TJ ' bl ( r ) = bl 'i rJ ■ 

j=0 j=0 



a i( T ) = 2 [M*> r ) + h{~h r )} , 



From flODD , (IB^j) and (E23J, we have 



t n (x,iV+l) 
(5.26) 



2m r(n + l)r(iV - ra + 1) 



x 



+00 /-(o+) 



a (r) + 



bn(r) 



u 



1 + -y) 9o(u,r) 



Nm(u--) 

3 \ u J 



r n e~ NT dudr 



u 



;-l) n+1 r(n + N + 2)e n b~ n e N ^ 
T(N -n + l)N n +^ 



c J (2Nm) + d Ji(2Nm) + 
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where 



(5.27) 



and 



iV" +1 r , s n - Nta T{n + m+l) 



r(n + i)y u ^ ' u "' l r(n + i)iv 



m=0 

oo 



JV"+i r „ ^ r(n + m + l) 



r(n + l)N r ' 

m=0 K 1 



(5.28) A = f°° / <0+) U + 1) ^il) e -(-ax»e--d M d T . 

2r(n + 1)tt2 Jo J_ oc V u J u 

In (15.261) . we have made use of the integral representations of Bessel function and Gamma 
function 

Uz) = hj^ ^" lexp {l( M -^)} dM 

and 

/"+oo 

m+1 / -Wt n 



r(n + 1) = iV™ +i / e- JVr r n dr; 

see [3 (10.9.19) and (5.9.1)]. Using (15.231) and integration by parts, we can rewrite in 
fl5J2HD as 

(5.29) 4 = / + " / (0+) 2^l) e -(»-a Tne - dudT . 

1 2r(n + l)7rii„ « 

Repeating the procedure above, we obtain 

(-l) n+1 r(n + N + 2)e n b~ n e N ~< r ' 



t n (x,N+l) 



(5.30) 



where 



T(N-n + l)N n +i 



Jo(2iVm) 



_Q_ 

N l 



1=0 

p-i 



-J l( 2iVm)^A + < 



z=o 



^ n+1 r-r-^-Nr^. ^ n I> + m + l) 



]\jn+L roc w 



m=0 

oo 



T(n + l)JV n 



(5.32) 



^ n+1 r h (\n-Nr A \- h r(n + m + l) 
N n+1 ~ p f + °° h p (u,r) Nm(u _ h) _ 



(5.33) £+ = — ^ — / / ' ?r,v e lv T-^" e-^d udr. 

2T(n + l)mJ J _^ u 



Since ho(u, r) involves df/du and dv/du depends on the parameters a and b in (11.31) . the 
coefficients C\ and d\ in (15.311) and (15.321) also depend on a and b. For an estimate on these 
coefficients, see (15.571) below. 
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To facilitate the application of expansion (I5.30p . we recall that the constants 7 and 
m are explicitly given in (I5.17P and (I5.18p . and note that the leading coefficient can be 
(asymptotically) calculated by using (I5.24p and (I5.25p . Indeed, we have 



c ~ ao,o = 2 i h o(i, 0) + h{-i, 0)] 

and 

do ~ &o,o = 2 M*> °) ~ K~h 0)] 
as b — > and A^ — ?■ 00. Furthermore, equation (I5.39P gives 



(5.34) fc o(±i ,0) = ^-^ J , 2{1 ~ 2a)m 

W oV4a-4a 2 -6 2 (l- % ) 



where A = 2 arctan ~~ VT^' 



5.1 The mapping v — > u in (15.151) 



In the case under discussion, a/b 2 — » 00 and x — > 00. Thus, since < a < |, we have 
4a — 4a 2 — 6 2 > 0, and the saddle points in (15.131) are complex. 

Theorem 1. When a G (0, |] ; b — )• + and 6 2 /a — >■ ; the mapping v — > u defined in 
h5.15\) is one-to-one and analytic for v G D v in the v-plane and u G D u in the u-plane, 
where the image of the boundary of D u in the Z -plane is given in Ii5. 36]) . and D v is the 
image of D u under the mapping Ii5.15\) . 

Proof. As in the cases of Charlier polynomials [9 J and Meixner polynomials [3], we intro- 
duce an intermediate variable Z defined by 

61n(l - t {v))+{l - b) lntoW - aln(l - v) - 61n(— -) - 7 

to{v) +v — 1 

(5.35) ( r 

=Z = m I u — 

where to( u ) is the relevant saddle point given in (15. 3p . 

We consider the upper half of the v-plane. Since the functions f(t (v),v) — 7 and 
g(u) are both symmetric with respect to the real line, the case for the lower half of the 
f-plane can be handled in the same manner. To avoid multi-valuedness around the saddle 
point, we divide the upper half of the v -plane into two parts by using the steepest descent 
path through f_, i.e., the point denoted by D; see Figure HI Call the two parts region 
I and region III. Note that there are two branch cuts: one along the infinite interval 
[1, 00), and the other along the line segment joining v = to v = 26 2 + i2by/T^¥, which 
was introduced in (15.31) . We denoted the point 26 2 + i2by/l — b 2 by the letter T in Figure 4. 

As v traverses along the boundary ABCDEFGA of region I, the image point Z tra- 
verses along the corresponding boundary of a region in the Z-plane; see Figure |5j In 
Figure [TJ we draw the boundary of the region, corresponding to region I, in the u-plane. 
The image point Z = ip{u) + 7 traverses along the boundary of the same region shown 
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Figure 4: The upper half of the v-plane. 







B 
C 


A 

Region I 

D 

G 




E 
F 




M 

N 


L 

P 

Region II 

S 


Q 

R 







Figure 5: The image of Region I in the Z-plane. 



in Figure [5] Note that we have used an arc GABC to avoid the cut OT . The boundary 
curves EF and GABC in the f-plane are rather arbitrary; for convenience, we choose 
them to be the ones whose images in the Z-plane are as indicated in Figure 5. 

From (I5.15p . it is readily seen that the mapping v — > u is the composite function of g~ x 
and f(t Q (v), v) — j. We have just verified that this mapping is one-to-one on the boundary 
of region I. By the same argument, one can prove that this mapping is also one-to-one 
on the boundary of region III. For the image of region III in the intermediate Z-plane, 
see Figure El In the w-plane, we let D u denote the union of the regions I to IV, with 
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Figure 6: The image of Region II in the Z-plane. 
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Figure 7: The upper half of the w-plane. 



outer boundary FEE'KJIVUTN'NM, and inner boundary GABCC HWQ'QRSL- see 
Figured Under the mapping Z — g(u), we express the images of parts of the outer and 
inner boundaries of D u in the Z-plane as follows: 



(5.36) 



GA 
AB 

BC and C'H 
IJ 
JK 
KF 



Re Z = bM 

-bM < Re Z <bM 
Re Z = -bM 
Re Z = aM 

- aM < Re Z < aM 
Re Z = -aM 



< Im Z < 6bn; 

Im Z = Obn; 
< Im Z < 9bn; 
< Im Z < an] 



Im Z 



air: 



< Im Z < air, 
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where M > is any fixed constant, and 9 is a generic symbol for a constant in (0, 1). 
Furthermore, we denote by D v the corresponding region of D u in the v-plane. 

By Theorem 1.2.2 of [13, P-12], the mapping v — > u is one-to-one in the interior of both 
regions I and III in the upper half of the f-plane. As explained earlier, the one-to-one 
property of this mapping in the lower half of the w-plane can be established by using 
the symmetry of the functions with respect to the real axis. Note that the only possible 
singular points of the mapping v — > u in D v are at v = v T . Since the images of these 
points in the w-plane (i.e., u = ±z) are bounded, the mapping is in fact one-to-one and 
analytic in D v . 

□ 



5.2 Analyticity of ho(u, r) 

We now investigate the function ho(u,r) = h(u,r) given in ( 15.21D . By (15.1 Oft and ( 15. 19ft . 
we have 

(5-37) ^ 

m(r - b)t{l - t)(t + v - l)(u - i)(u + i) [(1 - 2a)v + \/ Ab 2 - Ab 2 v + 



h (u,r) 



v 2 



2a(l - a)r(l + b)(t - q (v))(t - t (v))u(v - v+)(v - vJ) 
when t b and u ^ ±i, and 



m(u -i){u + i) \{l - 2a)v + y/Ab 2 - 4b 2 v + v 2 
h {u,T) = 



2a(l — a)u(v — v+)(v — v) 

(5.38) , — 

%{v){l-t Q (v)){t {v)+v-l) 



X 



bsjAb 2 - Ab 2 v + v 2 



when t = b and u ^ ±z, where u G D u and D u is defined in the previous subsection. 
Note that u = is excluded from D u . It is easy to see that (t — l)/r is bounded by a 
constant, independent of b and a, as r — > for v G D v and u G D u . Therefore, h (u, r) is 
analytic in r for r in the neighbourhood of [0, oo). Moreover, h (u, r) is also analytic in 
u for u G D u , since the only singularities u = ±z in D u are both removable; indeed, we 
have 



(5.39) h o{ u,r)- 0- W -*>(* + <* -1> / 



r(l + b)(t- 4(v T ))(t - t (v T )) y (1 - v T )bV4a - Aa? - b 2 
when r b and u = ±i, and by the equation following (I5.19P we also have 



(5.40) h (u,r) 



2mt T (l-t T )(t T + v T - V 



b 2 v T {l - v T )y/Aa - Aa 2 - b 2 



when t = b and u = ±i. To reach ( I5.40p . we have written u± as u± = \/u± and move u 2 ± 
inside the square root in dv/du at u = u±\ see (I5.19p . 

For the analysis to be used in the next subsection, we now give an estimate for ho(u, 0) 
when u lies on the boundaries of D u , i.e, the outer boundary FEE'KJIVUTN'NM and 
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the inner boundary GABCC HWQ'QRSL shown in Figure [7J For convenience, let us 
denote the outer boundary by Co and the inner boundary by Cj. 



First, we show that for u on the inner boundary Cj, 
(5.41) 



_ m m 
Ci— < u < C 2 — , 



where C\ < C 2 and C, C±, C 2 are used, here and thereafter, as generic symbols for 
constants independent of u, v, m, a and b. Recall from ( 15.18)) that m ~ frarctan a/yz^ a s 
b — > 0. To prove (I5.4ip . we take the part AB of the inner boundary Cj as an illustration. 
Using (15.361) . we obtain 



(5.42) 



m 1 



1 



Re u 



< bM and m 1 + 



1 



\u\ 



Im U = ^fe7T. 



Since |Im w| < from the equality in (I5.42p we have \u\ 2 — ^\u\ + 1 > 0. Therefore, 

9bir 



(5.43) 



u < 



2m 



9bn 
2m 



1 < 



Cm 
bn 



Rewriting (I5.42p gives 
m 



(5.44) 



\u\ 



Re u\ < bM + m I Re w and 



m 



\u\ 



llm m| < Obn + m llm u| 



which lead to -S- < 6M + #6^ + m|u|. Substituting (I5.43P into the right-hand side of 
the last inequality, and combining the resulting inequality with ( I5.43p . we obtain (I5.4ip 
immediately. Similarly, for u on the outer boundary Co, we have the estimates 



„ a a 
C x — < \u\ < C 2 - 



(5.45) 

m ' ' m 
Hence, it follows from flCTj) . and (E37J) that 

(5.46) \h(u,0)\<C for«6C ; UC . 



5.3 Error bounds for the remainder 

To prove the asymptotic nature of the expansion in (I5.30p . we need to give precise esti- 
mates for its coefficients q, di in f )5.3ip and f)5.32p . and the error term e+ given in (I5.33p . 
since the derivative dv/du in ( I5.19P may blow up as v± approach 0, just like the case 
in Section 4. Because the coefficients q and di are related to a p (r) and b p {r) by (I5.24p . 
(I5.3ip and (I5.32p . let us first estimate a p {r) and b p {r). To this end, we define recursively 



A {u) = — —, B (u) 



and 



if, i \ 1 r i d 



A '+ M = ™( 1 + ^J \u + TJ A * iu) < 

If. 1 V 1 f 1 d 
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for p — 0, 1, 2, .... By induction, it can be shown that 



(5.47) 



(2m) (it) (l + u 2 ) S d ^(i + u2 ) 



for k = 0,1, 2,..., where c^i and d^/ are constants independent of m and w (but dependent 
of r); see [12]. Furthermore, by (j5.25p . (I5.22p and (I5.23p . it can be proved that 

(5.48) a p {r) = ^- h {u,r)A p {u)du 

Z7rt JCnllCr 



and 



(5.49) b p (r) = ^- [ h (u,r)B p (u)du; 

see also [12]. Define q := min{a, b}; since in this case both n and x — > 00, and b 2 /a — > 0, 
we have qN — > 00, m/q = 0(b/ y/a) — > if a < b, and m/g ~ arctan \Ja/{l — a) if a > b. 

Using (15.471) . (I5.4ip and (I5.45p . it is easily verified that 
(5-50) \A p (u)\<^, \B p {u)\<^ 



for it on the inner boundary Cj, and 



Cm , „ , s , Cm 
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(5-51) |4, (tl )|<_ |S p (r*)|<^ 

for m on the outer boundary C , where C is used again as a generic symbol for constants 
independent of u, a and b. Also we have the estimates 

5.52 / ds < — — and / ds < — . 

Jcj b J Co m 

To see this, we take one part of Ci, namely, AB as an example; see (I5.35P and (I5.36p . 
where Z in (I5.35P satisfies Im Z = Obn. Using the second equality in (15.361) . and in polar 
coordinates, u = re tlB , the last equation is equivalent to 

(r 2 + l) „ dbn 

5.53 v ; sin/3 = . 

r m 

Thus 

(5.54) dr= \ - J cos/M/3. 

(1 — r z )6b'K 

Since in this case r < 1, we have |dr| < |d/3|. By (15.411) . |r| < C2fn/b and 

Cm 



(5.55) |da| = vW+^W < -y- |d/?| 

26 



from which the first inequality in f !5.52j) follows. The second inequality can be established 
in a similar manner. 



Since h o (u,0) ^ and h (u,r) is analytic in r, we have \h (u,r)\ < C\h (u,0)\ for 
some constant C > and for r in the neighbourhood of the origin. By a combination of 
Q53SD, d53S]), fl530|) . ([ESI]), f l532|) and fl5^6|) . we obtain 

(5.56) |op(r)|<^ and |6 p (r)| < j p 

for r near r = 0. Furthermore, by (I5.24p . (I5.3ip . (I5.32p and (I5.56p . we have 

C C 

(5.57) IcJ < — and IgLI < — . 
v J i pi - qP i pi - qP 

To estimate the remainder in (I5.33p . we split the loop contour into two parts; see Figure 
El The bounded part of the contour, denoted by r 1; is contained in a subdomain D u of 
D u , which has a distance c\a/m from the outer boundary of D u (i.e., Co), and has a 
distance c 2 m/b from the inner boundary of D u (i.e., Cj), C\ and c 2 being two constants 
independent of u, a and b. Note that c\ajm may become large, whereas c^vnjb may 
approach zero. The unbounded part of the contour, denoted by T2, is the rest of the loop 
outside the subdomain D u . Put 



(5.58) eU = ^ r [ M^ e M^) T v--d.dr, 
1 ^ pA 2T(n + l)mJ J Tl u 

(5.59) < 2 = ^t^- [ + °° [ ^^le Nm ( u -^T n e- N ^dudT. 



2T(n + ljiri J J T2 u 




Figure 8: The contour in w-plane. 



As in [5] and [12], it can be shown that el 2 is exponentially small in comparison with e 



27 



To estimate ei^, we also follow j5] and [12] by first deriving an integral representation for 
h k (u,r). Define recursively 

Qo(u,w) 



u — w 

(5.60) 



for p = 0, 1, 2, • • • . Using induction, one can easily verify that 

p-l p-i / 1 \ p+i / 1 N 

(5.6D Qs(B , m)= _gx: fej , j (_j „^( 



li — 10 

where g^jj are some constants depending only on i, j and p. Similar to (14.36p . for w G .D w 
we have from ( 15. 22ft . (15. 23ft and ( I5.60P the integral representation 

(5.62) h p (w,T) = — I h (u,T)Q p (u,w)du. 



CiUCq 



From (I5.6ip . it is easy to see that 



\Q p (u,w)\ < C — — for u G Co and w G D u , 

(5.63) a \ 

\Qp(u,w)\ < C — — — for u G Cj and w G D u . 
mb p 1 

Coupling the above inequalities with (I5.52j) and (I5.62p . we obtain 

(5.64) IV^t)! < ^ 

for w G D u and r in the neighbourhood of r = 0. 

If Nm is bounded, one can easily show from (I5.58P and (I5.64p that 

< 5 «> Kil = {^Ni 

If Nm — > oo we divide the integration path r x into three pieces: the steepest descent 
path L\ through u = u + = +i in the upper half of the u-plane, the steepest descent path 
L 2 through u — tt_ — —i in the lower half of the w-plane, and the circular arc \u\ = 1 — 8q, 
denoted by L 3 , which joins the steepest descent paths Li and L 2 at u = and u — uj, 
respectively; see Figure El Hence, 

(5.66) e+ A = h + I 2 + h, 

where Jj, i = 1,2,3, denotes the integral over the subcontour Lj. Applying the steepest 
descent method [TH p. 84], and using ( I5.64p . it can be easily verified that 



(5.67) |/i, 2 | = O 



1 



qPNP V Nm 
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For I3, since 



1\ (2 — <5)<y„ (2 , 



Re ( u 

u 



< Re I u - - 



,± \ u 



u=u± 



I2 is exponentially small in comparison with Therefore, 

1 



(5.68) Ki\=0 



p ' \qPNP\INm 

It is well-known that the Bessel function J a (t) is bounded when t is bounded, and that 
as t — > 00, 

(5.69) J a (t) ~ \f^~ t cos (t ~ \ a7( _ ^ 7r 

Coupling the estimates (I5.65P and ( I5.68P , and together with (I5.69P and (I5.57P , we can find 
a constant M independent of a and b such that 

(5.70) \e+\ \J (2Nm)\ l -^4 + \ J^Nm)] ^ 



which establishes the asymptotic nature of the expansion in (I5.30p . given that qN — > 00. 
In conclusion, in the case n 2 /xN — > and x — > 00, f)5.30p is an asymptotic expansion. 

6 REMAINING CASES 

We are now left with only two easy cases to consider. 

6.1 a > 0, a/b 2 large and x bounded 

In this case, the series in (11. 2p is itself an asymptotic expansion, since 

MN):= (-N + l) k k\k\ ' fc = 0,l,--",n, 

is an asymptotic sequence when xN/n 2 — > 00; see Table 1 and also [HI p. 10]. 

6.2 a < 

The case a < is relatively easy, compared with the case a > 0. Let us start with the 
integral representation (12.231) . For fixed w 6 72, the phase function f(t,w) in ( I2.23P has 
a saddle point 

2w-l+ ^l + Ab 2 (w-l)w 
(0) to{w) = 2jl + bj^ ; 
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cf.f l23|) . Note that the only difference between f(t, w) in f[2T2]) and f(t, w) in f[2T24"]) is the 
choice of branch cuts in the iy-plane due to the term alnw — a\n(w — 1) in f(t, w) and 
the term aln(— w) — aln(l — w) in f(t, w). Similar to (13.11) . for any fixed w G 72, we have 
to(w) = 1 — b + 0{b 2 ) as b — »■ 0; in particular, to(0) = 1 — 6. By the same reasoning given 
for (13.21) . we introduce the mapping t — > r defined by 

(6.2) b\nT-T + A = f{t,w)-f{t {w),w) 
with 

(6.3) r(t H) = b. 

Coupling (16.31) and (16. 2p yields A = b — 6 In 6. From (I6.2p . we have 
dt b/r-1 (r - b)t(l - t) [1 - (1 - t)w] 



dr df(t,w)/dt T(l + b)w(t-t+(w))(t-t (w))' 

(6 ' 4) dt = f t + H(i-t + H)[i-(i-t+HH ' 1/2 

dr [ 60 + 4b 2 (w - l)w 

where we have used L'Hospital's rule for r = b and 

2w - 1 ± y/1 + 46 2 (w - l)w 



r^, 



(6.5) t±( W ) 



2(1 + b)w 



cf. (l3.5p . Furthermore, it follows from (I6.2p and (I2.23P that 

(-1)™ T(n + N + 2)e n b- n 

UX ' iV + 1) - 2m T(n + mN-n + l) 
( 6 - 6 ) r+00 , 1 . df 

72 w — 1 dr 
cf. (l3.6p . Solving the equation df(t (w),w)/dw = 0, we obtain the saddle points 



b ± Vb 2 - 4a + 4a 2 
w± = 



(6.7) 2?; 



1 1 / 4a(l - a) 



2 2V 6 2 

cf. (l2.9p . The relevant saddle point on the integral path 72 is the negative saddle point to_. 

The following procedure is the same as that given in [8]. Recall the Hankel integral for 
the Gamma function 

1 1 

(6.8) = e u u- z du, 

T(z) 2m 

where the contour is a loop starting at —00, encircling the origin in the counterclockwise 
direction and returning to —00. With z = — aiV + 1 and u replaced by —u, we obtain 

(6.9) = — / -e N ^du, 

K ' T{-aN + l) 2m u 
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where ip{u) = aln(— u) — u. Make the transformation w — >• u defined by 

(6.10) f(t (w),w) = a\n(-u)-u + 7 
with 

(6.11) u{wJ) = a, 

where 7 is a constant to be determined. We have from (I6.10p 

dw _ (a - u)w{w - 1) [(1 - 2a) + y/4b 2 w 2 - 4b 2 w + 1 

(6.12) 



du 



2ub 2 {w — w + ){w — wS) 

dw f (l-o)(l-2a) 1 1/2 
du \ bWb 2 - 4a + 4a 2 J 
where we again have used L'Hdspital's rule for u = a. 

By (gElD and dUinD , we have 

/ ,r \ (-l)"T(n + iV + 2)e n b- n e N ^ 
t n (x,N + l) = y 

(6.13) 



u 7^ a, 



x 



2vri r(n + l)r(JV - n + 1) 

+ °° A ) r) jv(aln(-u)- u 



Ar(aln(- U )-«) e Ar(61nr-r) dudT 



where 
(6.14) 

Define recursively 
(6.15) 

(6.16) 



J 00 



h(u, t) 



u 



u dt dw 
w — 1 dr du 



ht (u,t) = ai (r) + (u- a)g t (u,t), 
dgi(u,r) 



hi+i(u,r) = u- 



du 



where }iq(u,t) = h(u,r), and expand a/(r) into a Maclaurin series 



(6.17) 



3=0 



T 



By an integration-by-parts procedure, we obtain 

(-l) n r(n + iV + 2)e n b~ n e N ^ 



(6.18) 

where 
(6.19) 

and 
(6.20) 



t n (x,N+l) 



T(N -n + T)T{-aN + l)N aN + n+1 



p-i 



Jfi + e P 



1=0 



ci 



T(n + 1) 



^° m=0 



r(n + m + 1) 
1 r(n + l)iV m 



iV a7V + n + l- Pr (_ aiV+1 ) ,+00 ,(o+) /^( Uj7 .) 



2r(n + 1)717 



./-oo 



3 ^(«) r » e -^dudr. 



Estimation of the error term given in (I6.20p is exactly the same as what we have done in 
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